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Pflm po tecnico 

La present invenzione riguarda m metodo per Telaboradone dei segoali 
. ecografici, iB specie nel campo degU esanu clirnd nor. invasm 
ultrasonica.L'kvenzioneriguardaanchemdispositivoecogc^^ 
10 im metodo dielaborazione da segnaHecografici. 
RtatftHftnatecnica 

,^a»^»<„.vive,rttrfteg»nche.Coine,«*>^ 

15 diu»as=riedi«g^ult««»ici(i»n„.»aunfi^^ 

.flitetenK, dd co<po dd pazie^ e ch. ecdtano m segaal. ecografico 

aetenxunato dalla riflesdone e/o dalla diffi>*»» aU'i^fi^ 

d« diffli««i e'o del riflettori Che « trowo aU-mt«». dcUa ZOB^ 

aieodsdoneoricezloBedelse^ecograficoix^sceadmsi^diricczi™^ 

20 elabo^ziooe U ^gode di ri^ die viene poi traaftm^ 

ecografica bidimensioiiale. 

GB ecogrf commetcidl foimscoBO i=»nagml relative aUWe^ ^ 
ecografico r*op«WP^ dai diverai tessud o interftcce investtte dal 

aciistico. 

25 I diversi tessuti o le diverse interfacce incontrate sono rilevate sulla base della 
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loro ecogemdta deterndnata dalla loro inq,edenza acustica che a sua volta 6 il 
prodotto Ira la deasita e la velocity acustica. 

Lo scqpo del tradmonaH ecografi h quelle di ricostruiie e mantenere 
rordinamaito spaziale del tessuto investigato attraverso an^dezze rivelate nel 
5 donmuodeltempocheasuavoltafoiniscelalocalizzazioneq)^^ 

coBdizionainento del segnale che viene e^ttuato. ha come obiettivi queUi di 
^na^^inuzzare 11 rapporto segnale rumore e di xnanteaere costante la sensibility 
dell'ind^iiie al variare deUa profonditL 

Quando ronda acustica si propaga in un mezzo ed incontra una zona 
10 disomogenea parte deU'energia viene xetropropagata verso il Wuttore e parte si 

propaga aU'intemo del tessuto. 

n segnale ,acustico retropropagato 6 costituito dalla sovr^sizione 
temporale dei conlribud provenienti da ciascun riflettore e/o difEusore; quello che ne 
risulta h un segnile con una forma estremainente complessa legata alia natura casuale 
15 della distribuzione, della forma e delle propriety meccaniche dei difiEusori e dd 
liflettori n segnale acustico retropropagato dal mezzo in esame viene convertito dal 
trasduttore in un segnale elettrico, denominate segnale a radio-ftequenza. 

In Fig.l h ripprtato un tratto di segnale a radio-firequenza appartenenle ad una 
tracciato ecografico relativo a tessuto biologico. m in parlicolare, la Fig.l mostra 
20 una porzionedi una traccia del segnale ecografico diritomo da miaporzione^ 

organo sottoposto ad esame ecografico. hx asdssa ^ riportato H tempo di risposta, 
cio^ Ja profonditi a cui si trova il diffosore od il riflettore che il segnale di 

ritomo captato dal trasduttore. In ordinata h riportata I'ampiezza del segnale. SuUa 
base di un in^ieme di tracce di questo tipo, che formano un fi:ame di segnali, si 
26 ricostruisce I'immagine ecografica dell'organo sotto esame. 
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Una prima elaborazione molto sempUce del segnale a radioftequenza 6 
I'estrazione deU'ampiezza del segnale mediante un'operazione di estrazione di 
invil^ppo. Nella Fig.2 h indicato con I il risultalo dell'operazione di invilnppo sul 
segnale a radiofrequenza riportato in Fig. 1. Si noti come il segaale di invilnppo 
5 segua le creste positive del segaale a radio-frequenza e pertanto dia mm stima 
dell'ampiezza del segnale. 

L'operazione.di inviluppo viene coirentemeate nsata in tatte le tecniche di 
imaging ecografico e sta alia base della preseotazione A-mode (dove A sta per 
amplitude) del singolo tracciato e B-mode (B dall'inglese brightness doh Imninositi) 
10 perleimmaginiottenutecodificandoinscaladigrigiodiboloier^^ 

singolo tracciato A-mode. 

NegU ecografi deUa attuale generazione si usa rinvUi^^po del segnale. 
ecografico, owero si r^presenta nei tracciati A-mode, B-mode, M-mode, fl segnale 
rivelato. 

15 Con questa soluzione si perdono le informazioni contemite nel segnale a 

radioftequenza, ciofe nella fomia generata dal trasduttore ricevente consistente in mi 
segnale con ftequenza portante di valore attomo a qnello della frequenza di 
. cenliobanda del trasdottore ricevente stesso, con mm modulazione di freqneoza 
(spettro asimmetrico) con salti di fese (picchi nella fiequenza) ed eventaalmente con 

20 morfologia.da distorsione non lineare quando il mezzo presenti effetti non lineari 

(nasdta di spettri armonici). 

In Fig.3 e riportata, in mia classica rappresentazione A-mode, il segnale 
invil^po del tracciato in radiofiequenza di Fig.l. Si noti come informazioni di 
ftequenza e fase del segnale originario siano andate iirimedi^ilmente perdute. 
25 Inparticolare, se si elimdna la radiofiequenza con mia operazione di estrazione di 
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inviluppo si perdono: 

. leMonnazioniderivantidaiBteraziomtmidtra^ 

non lineare (impedenza acustica locale ohe dipende daUa d«as^^ 
sua volta viene. a dipendere dalla intensit^t dell'onda ^Irasonica incidente, 
5 liemissione sulle armomche e sub armoniche da parte delle bolHdne che sono nel 

mezzo di coutrastoa ...etc) 
. ieiirfonnaziomderivantidaiBteraziomditipolmeare,cioHem 

spettrb del segoale per causa dell'assorbimento selettivo soprattutto in 
dipendenza dalla prpfonditS, I'effetto di scalloping (firastagUatoa della forma 
10 dello spettro di ampiezza presente in segnaK cha sobo' stati retrodif&isi da 
discontinuity anche puntifoniu distribuite periodicamente neUo spazio) 

I sistenu ecografici convenzionaH utilizzati per diagnostica medica si b^^^ 

qdndi, sull'osservazione di m tracciato A-mode o di una immagine B-mode in cui 
sono messe in evidenza esclusivamente le ampiezze del segnale eco relative alle 
15 diverse strutturetissutaliincontrate. 

Esistono,dunque,situaziomincailatecnica tradizionale mosim i suoi limiti 

in quanta, o per scarsa sensibiUt^^ o specificity, non h possibile giungere ad una 
diagnosi certa di patologia. Un segnale a radiofi^quenza, invece, contiene non solo le 
infonna^oni relative all'ampiezza, ma anche quelle di ftequenza e di fese Gegata ai 
20 ritarditemporaHconcuigliecMdovutiasingoUdiffusoriraggiungonoiltirasdu^^^^^^ 

Le ultime tecoiche di elaborazione del segnale tendono a estraire tutti i 
parameiri che caxatterizzano il segnale, e in particolare il cosiddetto contenuto 
spettrale. L'analisi spettrale appUcata direttameote sul segnale a radioftequenza^ 
' grado di estrarre iuformazioni da tutte le caratteiisticbe del segnale eco cosi |g^e 
25 essoritomadaltrasduttoie. 
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^ ddta penned di disport, di valori camtterisfici utOi alia 

differer^one ddle diverse tipologie del segnale e delle strutture 

biologiche con cvii interagiscoiio. 
5 B,Kg.4neldiag.anm».Mc.tocoa(A)4.iH»tatom^va.neateilsegn=lea 

«dioire^diFig.l.mertreidiagra.»ri(B)e(Qmo,.xauo.perd^ 

di Fig.4A, il oorrtea* speSxale del ae9>ale, p« evide^dare come le 
diversitt sWturaU nei tesauli porfeo a difiereoti tipologie di in«»bo eoografioo e 

.^.....^.....^^ ad una pe»^ 
10 spearo-ncor^o^tospettrale dene doe po.do»idi.eS«ae-«w^ 

4C si esteode da 0 a 12 MHz. Si Che enwn*i gU spetrt poadedono aU-i.^ la 

^ ftequeoza eeotrale (ftequenza in eorrisponda^ de^ 

^•inlonro di 4 MH. P seg^ale . reladvo ad 1^ scansione con sorrda .^.^ 
^ fteque,^ Bontale di 5 MHz). ™ preserte^, «n diverse a«ianiea» 

15 fierpeaziale,iB,u<^nel prime grafleoso.^1^ 

d« so«. pres«m irrr-toramil*. Ktore .KB'altra porziorie di segnale. 

In pratica, i seg-ali ultrasenici rellodi£M dalle nucrostr«iu.e, per le 
diffi^ di percorso di propagazio.«, i*ferisco.o aulla ftccia del »sd*or. 
ricev«* in n^adera distt«iva o ««.rumva da.»io luogo neUe inm^gW al 
20 fe.»me«oeosiddettodi^..L-ampiezzadeisegaaleeoo.do*reeo^^ 

eleva^ a seoonda del «po di interfered * anehe influeazaia dalla rifle«ivi« 
deitessudcheasuavoBa^de.— dalla«at„rameccanie..demiaedelasticia, 

aa eon^ lamiorostru^^ ad esempio un alto eon.eau.0 di ooUagene amenta 
rimpede.^ meccanica locale e ^ I'eoogemaa Per.a»^ aU'aspe*. 
25 deU-in=n»ginc»l«so»caclaasicaconM,»scopo:l-agglon^ones^^ 
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Che le e la !<». 

cha I'-gedda L-a.d™>«..o ddl'an^e^a media ndlo 

^„^tad».,»e»nfbn^«Baoa««edzz,^<-«nh,ocadei.^ 
U^u^o^rtadadiversifi^rimorfoMdchesicompongoBoinn^riemale^ 
Se si rie»» ad este-rte infbnnaziom suUo ^.tte in «nnW di 

^ (fte^ spaziale) si ac^uisi^- tof"-^- 
periodida spaziak delk aggjomeiaziom sotto esamo. 

n segpalc ecografico a radiofe,<«oza \ infl^ daU-infi«.a^«» che 
^ ^c* dene imn^gim. i. «mi« di pa^icolari o^dulazi.™ 

di fase del segnale ecografico. 

Nel donuMO della ftequeaza A ooaservano le infcnnariom prea^fi "el 
a<,^„ del .eo-po . ri^ d^o si tradue. nen-.ffe«o scal,^. che 

^i.™avanazienediampiezzas=o<»»lolascheo^o»g«&ari^ 

Fig4D. l.«d. ad es^n^o ■el oaso di N diffl^ri o nfle«ori CU. disHx^^ 
,.al^nonuoivooa.e«ain«la.i»oealUi™ghezzad-<^delsegnaled^ 

dal ftscio al^asoaico la de.^ti s^e di pote^ caloola.a secoado la 
a Fo^e.. de, se^ eoo ricevuu. pu6 essere <a«»eae««a 



come: 



Y^r^ +2^ S'-zO cos^^fc 



3 

dove: 1^ : density spettrale di potenza del segnale trasmesso 
ri: coeffidentediriflessionedeli-eamoriflettore 

U: distanza tea U primo e 1' i-esimo riflettore. 
c: vdodlii degli iilteasuoni. 
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Nel caso di sou doe iffldtori separati da uBa diste^^ pari a U'espres^ 
tma forma pii semplice: 

r . i?L = 2/ ^ la distanza nel tempo tra i pic^H dei segoaU retrop^^ 



riflettori. 

Si pu6 vedere dalle fonnule scheinatiiche che lo spettro tisulta pesato per un 
coseno con periodo propomonale alia distanza dei riflettori; questo detemuna dei 
valori di minimo e di massimo nella density spettrale di potenza. Questo effetto t 
10 cluan^.c.//oi,/ng.Taleeffetto6rappresenta^ 

. riflettori, posti a distanzevariabm (indicate con So. SI, S2)l'imoris^ 

nella diiezione di propagazione del fescio nltrasonico. 

Per poter distinguere due riflettori allrimenti indistinguibili, cip6 a distanza 
n,iaoredeUadnratadeU'impulsotrasmesso,sipuovalutareladist^ 
15 nel donnniospettrale. Nel caso dipia riflettori 10 spettro del segnale^ 

risulta pesato da oscillazioni sovxapposte di vario periodo. Nelle formule compaiono 
anche i coefficient! di riflessione dei riflettori che determinano I'ampiez^ del 

segnale retrodifiEuso. 

Dalla formula si vede che le informazioni present! nel dominio del tempo 
20 - influenzano la forma deUo spettro, per cui un'accurata descrizione dell'andamento 
della distribuzione delle axnpiezze in questo dominio dovrebbe fomire gU elemenii 
per la differenziazione delle stratture biologiche. 

Da queste considera2ioni esemplificative risulta come I'anaUsi spettxale sul 
segnalearadio-frequenzasiaunpotentestnmentoperlacaratta^^ 

25 medianteilrelativo segnale ultrasonico. 
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NegU -iM »ai la me^ a pm.i. di .ruovi algorteB pet relabcnrazLon. M 
^ ha consentito di ril«va« Wonnariom ulteriori dei .ess«4 andando ad 
^ oc. a &»io :*a»mc<, to ten»a di variazioM spemaU 

todrtte sul segnale riflesso. E' state. vcrificato le sWmre paU.Iogid» 

selettivo in frequenza. 

Le metodologi. sin qui svitappate. tottavia. non ^ risaltate efficiend, vooi 

„on idmee a carattert^a- in »odo afflddite il tessoto^ 

pesantezzi co=*aU>zionale degU algori.^ uffiz^ - 
10 in^,«^,3^<machepermemUviauali.^inteS^5^deiri^ 

prey=d«e =wU<-d<rf dinich. dove U v.l<«iti di 
fondamentale. 

m EP-A-1341003 e descritto un nuovo metodo. per I'anaM spettrale di xm 
se^e ecografico, che supera i lixniti e gli inconvenienti deUe tecniche aBteriori e 
15 Che consente di estrane l'i^o™e Che g^era lo 

conten^to di questa preceded pubbUcazione . inte^ente incoxportato nella 
presente descrizione e ne forma parte integrante. 
^ Rntnmflrio d^^lTTnvenzione 
Parteado dal metodo descritto in EP-A-1.341.003, scopo deUa preseute 
20 invenzione^queUodiformreunperfedonamentoedunampHamento^^ 

offerte da. questo tipo di analisi spettrale, per otteaere maggiori e v^M^ 

infermaaonisiillastnitliirasottoesame. 

fc««toza,partendodaUestessept«nesse alia base del metodo desS^;^ 

■ EP-A.1.341.003, mveodone qui descriSa oonse-rte di afBnar. la descri^ 6^^ ,^^^ 
25 spe«z» «^ >™ ''PP"*^ 
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o^W^^V'^'^^ inn>ododaestrairedal<:ampionesottoiadagine 
a mmero ffl tofcrmazlani possMU e/o ma pa effioaoe uHHzzazion, dri 

coeffidead <*e sono caratteiistici nella descrizione deUo ,>ettn> locale. 

to p«aca, s«»»<Jo primo a^. I'ta™^ * basata suU'idea di 
gU =*imalod locali. qui di segmto dmm^ 
locali, ricavaMi dal segnalc eoografico tramite »a procedirD««o analogo a qaello 
aesoiuo in EP-A-1.341.003, utilfazarfo mi-infcm>arioae relattva alia fcnna degU 
di distrib^zione in van. classi di valori degU *aaU,ri locali stessL la 
prafica, da ua flame di segoaU in radloftcq^eaza raccolto d*™^ struts 
si ricavano «ao o pi4 ftan« o n^d di paramcM ^ettaU locaH o stb^ 

Successivameate. si «ep» »na daborazion. sWisfica di ,«esti stimatod, 6c»ndo 

«^ «^ finestta di dimensiom p^deteoninate hn^ U n»Wc^^ 

loodi per d*=muna« ristogramau. a dismlmioae degU 

varie classi * valori. Da qu«ti is«g«aum si ricava nn paiamelB, legato alia lo» 
fonna. ad es=n,*o la deviadone standard. GU stin^tori looaU vengono poi 
combinaa. ad esempio moWpBcati, divisi od in ateo modo daborati insieme ai 
di «»ma ddl-istog«mma di distribuzione. « d>e si ottiene * ana nu^ 
^ee di sdmatori IccaU pesati. d» cons«^ di assodare inten«m di ™lori degU 

.amatori stessi a zone omogenee del campione invesdgato. TaH zone samnno poi 
individnate dali-aaalisi istologie. la mppres^-tazione di taU zone dd feme 
originale a^ene per mezzo di an codice eromadco so«appos.o Ha 

tradizionde iimnagine B-mode. 
' n metodo della presente invendone p«A esaere denominate sinteticameate 

i nietodo RUm: -IRadiofteqaency Ultrasomc 
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sull'impiego di stmmtori lo<^'che contengoxxo infonnazioni spettraH sul segnale 
idtrasonico in radiofreqaenza resdtuito daUa struttura esaminata. 

Come appariia chiaro daU'esposmone che segue, gli stimatoii locaU sono 
r^preseutati dai coeffideati di m polinomio interpolaxite di grado oppottuno, ad 
5 esen^iodiquartx,oquin1o grado. Che approssiir^i^^^^^ 

bande di ficequenza in cui 6 stato precedentemente suddiviso il segDale in 
xadioftequenza. digitalizzato e caxnpionalo. In altemativa. gU stimatori locali possono 
essere cosdtuili da combinazioBi dd coefficienti di vario grado del polinonno 
interpolante. La combinazione dei coefficienti di yario grado pu5 awenire prinm o 
10 dopo aver pesato dasoun coeffideate con il fettore'di fotm^ dell'istograxmna di 
distribuzione. 

Secondo un ^verso aspetto, I'invenzione comiste nell'elaborare i coefficienti 
del polinonno interpolante per ricavare stimatori IdcaU da una combinazione di due o 
pii coeffidenti di grado diverso dd polinomio stesso, senza pesare i coefficienti per 
15 il parameiro di finma deU'istogramma di distaWone. Anche in questo caso si 
ottiene, rispetto ai metodi noti. un maggiore conteauto informativo. in quanto si 
sfruttano piii coefficienti dd polinomio ixiterpolante tm loro combinati, anzich6 un 
coefficiente unico. 

La suddivisione in bande. do6 la decomposizione in bande o sottobande pu5 
20 awenire con criteri vari. Inoltre. andie il parametro di forma deU'istogramma di 
distribuzione degU stimatori locali pu5 essere scelto in vario modo. In particolare, si 
puo prevedere di 

1. utilizz^ diverse Wavelet Madre per adattarsi alio speckle tipico del tessuto 

I 

indagato; 

25 2. usare diversi UvelH di decomposizione sia per avere diverse risoluzioni 
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• .paziali, sia per utiH^ in maxdera combinata diverse suddivisioai deUo 
spettro complessivo a disposiaone. 
3. utilizzare comWBazioBi dd coefficienti del polinomio interpolate '-pesati" 
. con parameM legati alia mtura deUa lore distribuzione stadslica tramite mo 
studio orientato degU istogranmd ateaverso paramefii ^^^^ 

- classipiiipopolate 

- deviazione standard 

- curtosi 

- mom.entidigradosiq?eriore 

- indice di sinunetria 

La combina^one dei coefficienti opportunamente pesati geaera, per ogni 
finestra (impostabile in dimensiomO dell'innnagine ecografica elaborata, ma nuova 
notice costituita da taU combinazioni. Su di essa si esegoe ma nnova statistica e si 
scelgono le classi associate alle differenti strulture cbe si desiderano caratteri^e. 
Per amientare la capacity di differenzia^one dei tessuti cio. I'individnazione delle 
classi "tipiche" si utilizzanuovamentem'anaUsicombinata deUe classi. 

Forme di realizzazione particolarmente vantaggiose del metodo secondo 
I'invenzione sono indicate neUe allegate rivendicarioxd e verramo descritte in 
„«ggiore dettagUo nel seguito con rif^^ento ad alcmi esempi di 

sperimeatali eseguite. 

Secondo »n divert asprtto, a»zioM eefettuare «m operazione di 
c^nbinazione o pcsatura t« gU »d«^ locdi ed i panm«M di fco-a degU 
Uog^ di dtea»zio=c degU samatori stesri. rinvo-rione prevede di e»«»are 
^ <:o«bi>^one di due o pi4 staatori locaH lo>x> . di utiHzza« ^ 
come sSmatore locale oonfl*ato. Ad esempio, si pui, pieved«e di 
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ogpi maUcc d»i coefficieatt di vario grado ad polinonno tatmpolarte. 

Succesrivamerte. ooeffidend inposHom omologhe di queste matdd pos»«. =asere 

,,,,rfA«ti.»les«x.lioio<«fficia.tidit^ 
5 teIorosonm>adomolti,JicatLSi<^ene.ind.arifiva.in^™doumma^^ 

di ooefficienS cia^mo dd quali 4 co^taito da ™a ft-zionc degli stoa«>ri locaU 

c<mteraia iieUe inaHd piecedmtemeiite 

Rreve descr '^^^^ dei disegjii 

n tovato vetrimagUo con^ s=ga«KJo la descrizioue e I'umto diaegno. a 
10 q^aemostam^praacafcrmadiattuarione^mlWa^addl-n^ 

particolare: le 

Figg.l a4A^D.gi^ini«cecedenzadesmtte.mostraaoimtiacd^ 
ed ^ r^sentazioBe del contenuto spettxale in due por^om del tracciato stesso; 
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Fig.5 nxostra uno schema a blocchi semplificato di m'^aiecchiatura 
ecografica; la 

Fig.6A mostra m segDale ecografico di ecdtazione inviato dal trasduttore 

nell'organo da esaroinare; le 

Figg.6B a 6D mostxano rispettivamente U segnale KF di risposta, la Wavelet 
20 Da.becMe.l6umizzatapermtrareUsegadeRFdigitalizza^ 

DWPT otteauddaimadecomposizioneinottobande;ie ^^c^' 
Figg.7A.7I mostano schematicameate la sequenza di elaborazi^^ul 
segnaleecograficooperatedalmetodosecondolapresenteinvenzione;^ 

Figg.7J, 7K mostxano schematicamente forme di attuazione modificate^ 
25 proceduradielaboTazionesecondol'invenzione;le /^i-H 
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Kgg.8A e 8B mostaDo le i»magii eoografiche B-mode ddU catotide e 
Figg.9Aal7B;la 

Kgg5A-9C mostrano I'andamento di polinoffli iitfetpolanS di qurto 
5 ^inw«parScola«aspUoazi»noaU-«m.di™aca»«de,»lativiallapc^e 

a placca lijndica, del sangue e di calcificaaone; la 

del poli»cmio tau^lante di quarto grado «Jattvi dla pam<»« di pta^ 

della carotide; la 

10 Fig.n mo*a gli i^gnmum dl di«ibaziane dei dorfMeoti del poli>^o 

interpolante nel saogue; la 

Fig^l2 mosto gU istog^nuni di distribodone dei coefficieati del polinomio 

interpolante nella zona di calcificazione; le 

Figg.l3A-13C mostrano I'andam^ dei polinomi interpolaixti del quinto 
15 ordine otteimti per la stessa appUcazione di Figg.9A-9C; le 

Figg.l4, 15 e 16 mostrano, analogamente alle Figg.10, 11 e 12, gU 
i^togranmu di distribuzione dei van coefficienli del polinonHO interpolante di qninto 
gtado nelle varie zone di placca lipidica, sangue e calcificazione; le 

Figg.l7A e 17B mostrano gli istogranmd di distribuzione dd coefficienti del 
20 polinonuo interpolante del quarto ordiue appUcato aU'analisi ^ 

prostadca; le 

Fig.l8A moslxa mo schema a bloccH di una cascata di ffltri perla 
deconxposizione del segnale in radioftequenza digitalizzato e campionato in sei 
bande o sottobande tramite una wavelet discreta, in bande di ampiezza variabile, con 
25 vma suddivisione piii fittanella zona delle basse ftequenze; la 
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Kg.l8B mostm aiemaltoameBte te sd tende di decoa^o«e: ta 
Hg.l8C nostra Vml^ dd. polinomio inteqpotate a qua* ordtoc 
„tte«^«ttca»doladeoon^onediHgg.l8A.18Badmseg.«lecco^^ 

risposta di xma carotide; e le 
5 Figg.l9.20e21mo*anocia»m™a>nm»«ineecograficaB-modediuna 

p<^c»e di oa^tide inte«ssata da placca Bpidl- c ft«>meni di oaIdfloazio«. 
raod,^ del po^^o interpotan. e gU istogranm. di distibu^one dci 

,tacoaHpMca,sa.,SOa(vasolft«o) = calcm^om,ri.pe«vam=n^^ 

^ . *nrtt,im ddl'albero di deSomlKMMone ottetrato ttamite 

10 Fig.22 mostra ima struttura aeu aiocru r 

ma trasfcnmta DWI; e le 

Figg^3Ae23Bmos.raoog«fidespUo.tWdenapro<:ed«mdidassfficazio^ 

de^ stimatori locali. 

« S,Fig.5itappre«>*in™do=s«o^«»d«maacomdiag^ 

Hocdd di uaa appareodua^™ eco^- con 1 * indicate oa 

di soan^one. genericamen.. indicaia e™ 2. cl« emcae il -gnalc "ogra^ 

ecci-azionc . riceve il segnale eoo^oo di rispo^ L'uscha analogica 3 deHateata 
a » eonegataad^aschedadi ac-luisiMone 5.Da<pesU.viene gene^ 

20 segaalea«diofte<,aenzacan.pian^edigi,aazzaU,.c.»eo^ 

di fiW FIR integral. i»&ato gen^iean^nte c<« 7. n l^nco di ffl 

^costoi*. daun qualnnque disposHivo integral e<««»rciata=a.e diaponibae, 
adesen^Uoireui*. 43168 p^daUaHa^a. a ci»ui.e OF191 prodot.0 daUa 
^ ■ o,„nn. ii dro* GEC 16256 p»de«e daPlessey o il dxcai.o GC2011 prcdo«o da 
^ O^eMp. n fflro HR 7 * as»d«o ad nn i^ceaso^ di segnale digiMe (DSP) 
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Micato con 9. ohe fa^ com. ccatroUo^ cto gesSc. I'intera ™nt» di 

dabcradooe. D ndcrocoolxollorc 9 pu6 es^ =d escnpio m. DSP moddlo TMS 

32031 delta Texas (USA)odaltroe4iiivaleiite. 

Con 11 4 indicate .ma BenK>ria associata al ffltro 7 ed al cortionoie 9^ 

5 n>en.odall*oonegataadmc„nverfito»analogico/digiMel3edaglis«di 

della testa di scansione 2. 

L„ sdneoa qm descritto i di per a6 conosduto e non lichicde deacrizioni 
pa^lari. In sintesi. I'app^o opera come seg»e. H Wn«ore I eme«e nna 
«de di segndi ecografid di ecdtazione aventi, ad eaen^io. randan»nto iUnst^ 

10 rig.6A, all-inlemo dd corpo dove si trova I'orgam-da esannnare. ad eaen.I»o la 
pros^ il fegrto od altio. n tcssuto invesfc dalle onde ^trasoniche dd segoale 
«^ dal WuBore 1 .esdWscc m> segnale difB.so o riflesso Che ha ta *^ 

segoale di Fig-l (segnale a rodiofiwpenza). 

a segnale a ladioftequa^ ^Hcne prdevato dalla catena di ricezione 
15 dell'ecoBrafo inmiediatamente dope U blocco di comj>«.sazione del guadagno in 
fcnzione del t«-^ (TGC) e q,teU digHalizzato con fieqneoza di campionamento 
cwortuna, ad esempio 40 MHz e con 12 W. (od atea risotazione adegnate) di 
risotazione dalla scheda di acquisizione 5. In praSca, I'acqmsizione avviene per 
^oletocceflno acompWamento dinnqnadro (o fian») „ di una porzione di 

20 .padro. I can>pioni di segnale acqm^ P« oiascuna singola traccia vengono 
...em^ in un banco di menK«ia delta scbeda di acqoidzione to ^ 

sta.0 acquisito «n quadro completo oppnre una porzionc di qoadro e solo a questo 
ponte i daa n^morizzatt >^ inviati alta sd-eda di elaborazione com 
DSP 9 ed a ffltro m 7 Che esegoe le operazioni d» desmlte in dettaglio 

25 nelseg«to.secondoilmetodo delta presenteinvenzione. Durante npa^ 
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dati ^ »*eda di dabon^onc. la sdKda di aoqui^one 

di segnale delle tracce reladve al quadro suooesrivo. H tempo di 
dalK^one della soheda di dab<»zi<«» 12 4 »fflci«ta«.te oontenuto da 
a con>pl*n««» ddl'elaborazl<m. ento a «^ « 

5 aegaalcc«^o„*«l=tt'oalq«^-««^--«'««^""^'"^^ 

ia tempo reale. 

L. sd>cda di elab<«azione 12 precede ad effittuare le elab<«aziom del 

segaale, ed in pardcoto la traaS»«-a tempo-fl«q«.« tr^ 

le atae »peraMom d-e ve^™ descritte in maggiore dettaglio nel scguito, per 

10 otKnerelaoaiatteriraajioiieddlospettrodelsegDato.' 

U Figg. 7A.7E riaasumoao in ft«na grafioa le operazioid s«bite dal 

in ingre^ per ottepere la rappresen^zio^ della maMoe ccn*^ 
^ve aU-opportn^ combinazione degH sdma«<m locaH detenninad .«nnte B 
^ „gge«o della i-ozione. Per nna d«crizio„e de«agliam deUc 

15 operaaio-i cegoite neUe vade M a elab<,raziane in nn esen^o .pedfico. si 
rin«oda a ,u««o descritto in EP-A-U41.003. ove non dettaglta^ 

questa sede. 

X^laborazione del segpale oampionato e digitalizzato ao<pnsto 
5 prev«l. in prin» hogo rappHcazione di nna taaft^a^ tempo-i^^ 

20 codi&areaseg.>aIein«mobande..n»..anza.peroiasc«»«cd^ 

ftame «xr«o. canipionaio e digitelizzato, U sef^e viene »=ompos.o in ^ 
<fi bande o so«obande di 6^ mediant l-appHcaziane di una tr^f^ 
^^n^fieqoenza. ad ese«pio nna trasltaata Wavele.. ed in pardcolare ai^pio 
.nediante la deoomposizione Wavelet Pacfats disoreta E fiteaggio a™ene «a^^| 

25 bana,diffl«in«egr«i7.nnmnerodiso«,bandeobandeincdvienedecon^P|^^^ 
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Ndl-esposmo,^ d» segue si ftri riferimealo ad ™m proccdum di decodifioa to otto 
bande o sottoba^ d» ve^no indicate Ml «ga»o oo=» ba«da 0, basda 1 .... 

baoda 7. Pe«teo, come apparixi chiaro in seguito. <pe^ valo» non ^ vinoolante, ed 
5 an aapcttoinm-vativo del n«*odooggrttoddla present invenzioM 

j^oaenacapacitidie^anedanostessosegnale^^ograficoinnuaoile^ 

^ maggiMe di infonna^oni per la caiatterizzazione tis^ude operando naa 
aeo„aificainnnnume.odiWeoso«»l>andedidiversalargi»^(cioiaiampieza 

variabae da banda abanda) e in uumero anche diverso da otto. 

10 I coeffidena d« filtri posaono esaere defiiti utiJiizando naa Wavelet 

Danbecbies 16. Eaen*lifloa.iv.Me.*e. in Fig.6C 6 riportato Tandamcnto della 
wavelet Daubeohies 16 ed in Fig.6D Tandamento del coefficiend DWPT o«eanti 
.fiteando il segoale RF di Fig.6B digitalizzato e campionalo. Peratoo, il t^ di 
rtafonn^a tempo-ii»penza ntilizzata pu5 variare aocbe in tozione del tipo di 

15 stn^o.easutodaanalizzare.Sipo»onoadeaen^iousareWavel«dh«se 

Danbednea 16. T^oamente. per I'e^ne della carotide. ad eaempio. pu6 es«« nsata 

^ Wavelet Wortbogonal 3 .7. per la caratterizzazione deUe placche carotidee. come 

verri eapostonel seguito con riferimaito ad alcuni dati apeiimertall. 

n mmiero dd coetadentt del ffltii i in genere baaao, ad esempio per a filtro 

20 Daub«=bies 16 . di 32. e risulta quindi ag^ole impl«n=ntarU su un oixcui.0 
o^mne^iale. che pu6 realize normalmente filM con ™ numero limi.a«> di 
coefflcient al massimo 64 per ogai sezione. con dedmazione per 4 e se a filtro non 

J a fee lineare. mentre con le d„e sezioni in cascata si pu6 arrivare a 128. I 
coefficienti dd fUri Che realizzano ina trasfbrmata Wavelet sono defini^ 
25 ricorsiva . sono di per s* conosciuti. Ad esempio posaono essere ntBizzati i 
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' coeffirie«.d.fl»MmUS-A-6(W6098,flcm««.=m.t6 4in^^ 

descrizLone. 

Jn Hg.7A 6 sdKmatiamcrte .WP«s«>teto fi^ 
testa di sca^ion. ddl'eoografo. H fiame. nd segato d=KmtoU> «ft«ne di 

ses^e amdogico raw^se^tato dal d-<»da 

ta cui si r<»da -toso^m di ecoiWm. Dopo a can>pionam«r,o e la 

digi.aliz.azi.me ri <«i=n. (Hg.7B) m .tkuJ™ 6 toe 0^ 

dene quaH coatiene serie di canapiom di .eg-^e digiMizza«.. 
10 ftaniB vet* indicate come «&amedigitalizza!o». 

ApplicaBdo la ttatfonnata tempo-fieqoenm al ftame digitalizzale. cio* 
il segmle can^onato e digitalizz-o costiSnto daU-i^cme dei coeffloienti 
a fiam. digitalizzato (scl^matizzato in Fig.7B), si o.ti«» 
in bande o aottobande: in pratlca per ciasc„na t«=oia da 1 a n si haim^^ 
15 istante di ten^. cie6 per can.pi<»« l»ngo la di-ezione dei tenrpi (cio» deUe 
j^ofcnditi di prepa^one del segnale), un coeffidente in ciascona defle We 
..ppesentale scton»«icamente in rw7C ed indicate c<»ne ^ 

banda?. 

NeU'esempio sd.«natizzato le bande sono otto . tnlK di nguale 1^^ 

20 Tnttevia.comesopraaccennato.sidevecompr«.dere(ean«W 

con rifertaerto ad alcmd dati sperinrortahO Che a nmn«. di 1^^ 
^ene scon^ U segrmle p«6 essere s»,«iare. Con m maggiore nnmero di 
sottobande di decodifica dal segnale in radioftequenza si o«ene una maggiore mole 
di informazioni snUe carattcrisdche spettraU a esse e qdndi snlla s.™^ 

25 indagine. tootee. la suddivisione in sot.>bande pu6 anche essere otte^ta con 
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intervalU di ftequenza (cio^ ampiez« di banda) non costaati, cio6 in bande o 
sottobande di ampiezza variabile e con filtri di forma qualsiasi. La scelta del tipo di 
deconyposizioneseguirioriteasperimeritaUinW^^ / 

struttura da sottopone ad indagine ecografica. 
5 La matrice contenente i coefficieixti deUa trasfoimata tempo-ftequenza viene 

indicata nel seguito come «matrice della irasfoxmata»; i coefficienti di detta 
n^ce vengono indicali come «coeffidenti spettxaH». Essa 6 una matrice 
tridimeasiox3.1e, le cui dimensiom sono date da: numero di tracce. numero di 
campioni per traccia, nmnero di baxKie o sottobande in cui e stato suddiviso il segnale 

10 inradioftequenza. 

L'operarione successiva consiste neU'estrarre gU stimatori locaU del segnale 
per ogni pnnto della maHice dei coefficienti spettrali (e quindi per i vari pmti del 
ftame di ingresso da cui la matrice dei coefficienti spettraU 6 stata ottenuta con le 

operazioni sin qui illustrate). 

A tale soopo. lm«o dascom. singolaliaeciawogcmo presi, pa: cfascun Wante 
cio4 po daaom. p»».o lm«o la din.e»done dd tempo (cio4 deHa 
poftnditt) i coeffidmB delle varie bar* 0 - 7 e VMM detamtato poltomio 
fatepota* cte approsA^a randamento di querf coeffidentL Questa aper.zi<me i 
«j««floa«. nelle Figg. m 7E. m m parti^lare, in Fi^^ 
20 ooffidenli spettraU neBe vane sottobamle,ielattvi alia tracdaled ato 
««„5,moath«merte ad m Waale teo^e pk quaUasi. T»W 1 ^ 
v»e ba,>le 0-7 in podzione U e dJ. sono riportafl tra loro alliBeaii mpetttvameni^ 
primo e nel secondo diagnunma di Fig.7B. dove con PI, 4 indicate U polinomio 
itopolante. oio* tt polinonno ehe appros^aa in modo otSmale i coeffidend 
25 spettraBallineadallaprof«>di.iiksullatiacdaleconP1.4indi^ 
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iBtei^^olantexelativo ai coefficieati spettrali allineati aUaprofondMk sullatracdan. 

Uno o m dei coefficieirti di questo polixxomio costituiscmo mo stiinatore 
locale Che cansente, come verr^ chiarilx, in seguxto, di estrane le i^rfonnazioni 
desiderate suUo spettro del segnale in radiofrequenza." Tale parametro, infatti, h 
6 indicative deflavaria^onedeUecaratteristichespettraUd^^ 

fetto Che esso e per definizione legato alia propriety del polinonuo che approssinm 
randameato dei coeffideoti spettraH in ogni banda per ciascun pmxto deUa mairice 
deicoefficientispettraliequhKHindirettamentedeltoeiningres^ 
lo stiniatore spettrale locale estratto dal polinonuo interpolante PI fomisce una 
) 10 indica^onequahtativadellavariaWlit^deUe carats 

in tm dfitenninato pvmto del frame in ingresso. 

Poich6 il poUnomio interpolante viene calcolato in corrispondenza di ogni 
p^to del frame ecografico di ingresso, a risultato di questa operazione 

eiascuno dd coeffidenti del polinomio inte^^^ 
15 dall'insieme degli stimatori locali. Questa malxice, nd seguito indicata come 
«marice degU stimatori locali>> o «frame degli stimatori locah» ^ rappresentata 
ux Fig.7F per un smgolo coeffiderrte del polinomio interpolante. Qui i singoU valori 
> .deUostimatorelocalesonoindicaticona«,doveivariadal an.essendonilnnmero 

di tracce della matrice, e j varia da 1 a m, esseado m il nnmero di coefficienti per 

20 iracda. Si deve comprendere che a seconda della significativit^ dd singoH 
coeffidenti di vario ordine dd polinomio interpolante si pu5 ipotizzare di cr 

matrice di stimatori locaU per dascun coeffidente dd polinomio. Per m 
interpolante dd quarto ordine. quindi. si possono costruire fino a cinque 
stimatori locaU, una per ciascuno dei dnque coeffidenti dd polinomio. I 
25 numerodimalrid si pu6a>struiieutilizzandoun polinomio digrado^ 
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,„,eas»iii.e.«e «ti i co«ffid«>a dd poltaondo posso»> essere utiU ndta 
...rt^izzazkme spettrale dd dr^olo segoal. eoogmfloo. U maggio« o minoB> 
sigdficattvia di coefBdeafi vi««. dcteBBmata =5«rm«tatac,rte in fmmone 
del tipo di tes^ o di strutuua da indagaie:Nd segmto ™m«™ 

esempi non Iknitativi. 

Per il calcolo del valoxe dei coefficienii del polinomio pu6 essere 
opportuaameate utilizzato il metodo di approssixnazione ai « quadrati. Tale 
nxetodo permette di detemuBare, dato m numa^^ 

y^, cormpoadenti alle ascisse Xo, xi, x„...il polinomio di grade m < n-1: 

P(pc)^ao + aix+ + am3i^" 

cbe approssiini nel miglior modo possibile i dati tilevati. 

Considerandp la sonuna S dei quadrati degU scarti Ira i valori ^prossinmti e 

quelli osservati: 

5 = 2(ao + a,xj + + a^xj -yj)\ 

affinche ii polinomio approssimi nel migUor modo. possibUe i dad, bisogna scegliere 
i coefficienii in maniera tale cbe 5 sia minimo. Derivando S ri^etto ai coeffidenti e 

e V. =yjc*v, si ottiene il 
uguagliando H risultato a zero e posto =±,Xj ® ^* ^ 

sistema di equazioni che in forma matriciale pu6 essere scritto come: 

dove a 6 U vettore dd coeffidenti dd polinomio P(x). 

Si 6 osservato dxe per vari tipi di tessuto e per varie corulizioni patologiche i 
. valori dd coeffidenti del polinomio interpolante assumono statisticamente valori 
determinati. Definendo classi di valori. si t inoltre osservato che la distribuzione 
statistica, cio6 I'istogramma di distribuzione di ciascun coeffideote del polinomio 
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inteipolante, presenta ma fonna variabUe a seconda del tipo di tessato e delle 
eventuaH condizioni patologiche di esso. Secondo nn aspetto della piesente 
invenzione, il metodo h basato su questa constatazione e ntilizza i parametei di foma 
dell'istogranuna per estraire m maggior niameco di infonnazioni spetixali dal segnale 
5 ecografico. 

Per comprendere I'importanza e I'entit^ di questa variazione della forma 
. dell'istogramma in varie condizioni e per vaxi tessod si feccia riferimento aUe Figg. 
daSAalTB. 

NeU'anali^ della carotide per verificare la presenza e la consistenza di 
10 eventuaUplacchecaiotidee.adesempio,^.emersocfiead»ttandonnad^ 

sedici bancle con ima Wavelet biortogonale 3.7 

ed titilizzando un polinomio interpolante del quarto ordine. non solo il terzo ma 
anche gli altri coefficienti del polinomio presentano una distnWone significativa 
per la caratterizzazione deilessuti indagati. Si e anche rflevato che nna analisi de^ 
15 istogrammi di distribuzione evidenzia mia diversa morfologia degU istogrammi di 
distiibuzione a seconda del vaii tipi di tessuto esaminato. 

La Fig.8A mostm I'immagine ecografica B-mode della carotide iitilizzata per 
oitenere i risultati sperimentali qni di seguito discussi. Snll'immagine sono anche 
individuate tre finesixe cenirate sulla placca Upidica, sulla zona hT)era del vaso, doe 
20 sul sangue, e sulle calcificadoni. Su queste zone presentanti sfrutture tissutaU diverse 
tra loro sono stad calcolati i coeffidenti della maMce degli stimatori locah, 
dimostrando che tali stimatori hanno nna distribuzione statistica variabile a seconda 

del tipo di tessuto indagato. 

In Fig.9A, 9B e 9C sono rappresentati gH andamenti dei tre polinomi 
25 interpolanti del quarto ordine relativi aUa placca, al sangue ed alle calcificazioni, 
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rispettivamente. la aliri tennuu, s^do il metodo di decodifica ed elaborazione del 
segnale. in radiofrequenza digitalizzato sin qui desmtto, i polinonn del quarto ordine 
Che approssimano in modo ottimale i coefficienti spettraH (cioe i coefficienti della 
ioatrice della trasfonnata) hanno I'andamento di Fig.9A per le porzidni di innnagine 
relative ad una placca Upidica, I'andamento di Fig,9B per le poizioni di innnagine 
• relative alia zona Ubera del vaso, in cui scone il sangue, e 1' andamento di Fig.9C per 
le porzioni di innnagine relative alle zone interessate da calcificazione. 

Nelle Figg.lO. 11 e 12 sono mostrali gU istogramnn di distribuzione statistica, 
nelle varie classi di appartenenza, dei cinque coefficienti dei polinonn interpolanti 
) delle Figg.9A, 9B e 9C rispettivamente, otteauti attraverso fanalia di diverse zone di 
^ porzione di carotide. In ascisse sono indicate le classi di appartenenza ed in 
ordinate le frequenze della popolazione. 

Per ciascuno dei cinque coefficieoli considerati (indicati qui con AO. Al, A2. 
A3. A4) sono slati costaiiti tre istogrannni com^ondenti a tre distinte finestre 
5 posizionate rispettiyamente: in una zona di placca Upidica OFig-lO), in lu^ zona 
Ubera del vasp,cio5pienadi sangue. (Fig.ll)ed in una zona soggettaacalcificazione 

(Fig.l2). Si osserva dagU istogranmd di Fig.ll, che nella z«na Ubera del vaso. si ha 
una distribuzione molto stretta dei coefficienti di ogni grado. I relativi istogranuni di 
. distribuzione presentano una dispersione molto modesta. Viceversa, neUe zone 

20 interessate dalle placche Upidiche e dalle caldficazioni (Figg.lO e 12). in modo 
particolare U coefficiente del secondo ordine (A2) presenta una distribuzione 
st^stica sostanziaUnente differente rispetto a queUa che presenta I'omologo 
coefficiente nel tessuto sano (sangue, Fig.ll). ed inoltre sostanziahnente differente 
nei due casi di placca Upidica e caldficazione. In altri tennini, anaUzzando 

25 I'andamento degli istogrammi rappresentativi deUa distribuzione statistica del 
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■ coeffloiertedelsecoU.<.rftaesi*mgradodidisBng«^ 

i in pt««.za di ™o Bbeto (sangoe, Fig.1 1). plaoca UpMoa (Fig.lO) o caldficazione 

ffig.l2). 

Q„«to rigmfica cto U placca Upiffica e la calcfficazicme 

5 segaalcie^WP^alt^aziorispetWidh^r^^ 

vaWaztaoe qualltadva deU-istog«.m..a di dtateWcme dd co^oie.*! dd «condo 
orftoe dd polmo.mo in»«potot= Qhe l-andamento d=i ooefficienti 

spettrali. 

Analoghe coBsidera^iom possono essere svolte suU'andamento 
► 10 deUMstograxnnm relative al coeffidente del terzo Srdin^ 

distribuzione diversa nelle tre finesde considerate, cenlxate su placca Upidica, sangue 
. e caldfica^ione rispettivameate. In misuxa minore, ma ancora significativa, 
differisconp gU istograxmxu xelativi aUa distribuzione del coeffidente del primo 
ordine. 

15 Utilizzando una decodifica in m elevato numero di bande, sedici 

neU'esempio dei dad riportali in queste figure, si possono utilizzare polinomi 
. interpolantidigradosuperioreal quarto, ad esen^iounpolix.0^^ 
> cpiinto ordine. Da questo si ottengono risultati sperimentaU ancora piii significadvi 

relativi alia distribu^one statistica dei vari coeffidenti. In Fig.l3A, 13B e 13C sono 
20 riportatigUandamentiddpoUnonuinteqK,lanti del quinto ordine cte 

randamento dd coeffidenti spettraU nd .caso deUa carotide. ottenuti con ^ 
decodifica in sedid sottobande tramite Wavelet biortogonale 3 .7. ^ o 

Le Figg.l4, 15 e 16 mostrano gli istograimni rappresentativi d^ 
dislxibu^one statistica dd coeffidenti (A0-A5) dd vari ordini del p^onuo^ 
25 intetpolante per lestessefinestre indicate in Fig.8A,centratexispet^^ 
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placca Hpidica, sul satfgue e su una zona interessata da calcificazioni. 

Da m confronto tra le Figg.9A, 9B, 9G e 13A, 13B, 13C si osserva m 
diverso andamento del polinomio inteipolante sc^ttutto nella zona della placca 
Hpidica e della calcificarione. A questo corrisponde im diverso andamento 
5 deUMstogrannnadeicoeffidentiddpolinoiniinterpolantineUevari^ 
' differeoza di distribuzione statistica tra i vad tipi di tessuti 

riq«tto a qnanto risconlxato negli istogranm^ 

Risultati di analogo interesse sono stali ottenuti neU'esame ecografico della 
ghiandola prostatica. L'immagine ecografica B-mode utilizzaia per ottenere i risultati 
10 sperinxentaU di seguito esposti e discussi 6 riportata"in Fig.8B. In essa h indicata 
anche la finestra centrata sulla zona di tessuto tumoiale esaminato. 

Anche nel caso della prostata, utilizzando un polinomio interpolate del 
quarto or<fine ed una decodifica del segnale a radioftequenza di^taliz^ in sedici 
bande o sottobande si riscontra una distribuidone statistica dei coefficienti del primo, 
15 secondo, terzo e quarto ordine del polinomio interpolante Che si diff^^ 

tessuto sano e tessuto neoplastico anche per forma degii istogrammi di distribuzione. 
Le Figg.l7A e 17B mostrano gU istogrammi deUa ftequenza di popolarione nelle 
varie classi di appartenenza per i coefficienti al e a2 (Fig.l7A) e per i coefficienti aS 
e a4 (Fig.l7B) del primo, secondo, terzo e quarto ordine rispettivamente, sia per U 
20 tessuto tumorale Che per il tessuto sano, oltre Che nel caso di assenza di tess^ 
(rumore misurato in un campione di acqua). Si osserva, dal confironto degli 
istogrammi, che tutti i coefficienti fomiscono una mfoimazione utile per 
rindividuazione del tessuto tumorale, in quanto le forme degU istogrammi dei vari 
coefficienti cambiano tra tessuto sano e neoplastico per tutti e quattro i ci>effidenti. 
25 Da queste risultanze sperimentaU emerge da un lato la rilevanza (ai fini 



C:\Documents and SettingsNimchele 



Pag.26di44 

081003\act300903def(2).doc 



deU'anaUsi spettrale) di diversi coefficieati del polmomio. interpolate, e daU'altro 
mportanza del fettore di forma degli istograxnnd di distrib^one nelle vane classi 
di valori degli stimatori locali ai m dell'individuazione dei lessud tmnorali o 
coimmque di particolari caiatteristiche degU organi sotto esame. 

Pertanto, secoBdo un primo aspetto. U metodo della preseate invenzione 
prevede I'utilizzo dei fettori di fonna degU istograaxuni di distribmione di mo o piii 

d^ coeffideod del polinomio interpolaate per otterxere ma informazioue 

rappresentabile in sovr^sizione alia nonnale inunagine B-xnode dell'ecografo. 

Secondo un diverso aspetto, H metodo della presente inven^one si caratteri^ per 

rin.piego(incombinazionetraloro)dipriistimatoril5^^^^^ 

vario ordine dello stesso polinomio inteipolante. 

Una possibile modaHta di elaborazione viene descritta nel seguitx, con 
riferimento aHe Figg.7G-7I. Nell'esempio viene considerato ma singola mabice di 
stimatari locali, cioe i valori di un solo coefficiente del polinonuo interpolante. La 
^natrice degU stimatori locali considerata pu6 essere quella relativa ad nno qualsiasi 
dei coeffidenti del polinomio interpolante, che ^ risultato sigr^cativo nell'indagine 
del tipo di tessuto sotto esame. La procedura che verri descritta nel seguito p«5 
essere ripetuta per ciascnna delle matrici degli stimatori locali di interesse, ogni 
n^ce essendo ottenuta raccogUendo i valori che assume, nell'immagine ecografica, 
^ determinato coefficiente del polinomio interpolante. Si pu6 anche ipotizzare di 
partire da una malrice di stimatori locaU costituita da una combinazione di matrici 
«aative a coefficient! di ordine diverso dello stesso polinomio. Ad esempio, si pu6 
prevedere di generate una matrice di stimatori locaH che raccogUe i coefficienti del 
secondo ordine ed una diversa mairice di stimatori locaU che raccoglie i coefficienti 
del terzo ordine e poi da queste una matrice di stimatori locali in cui ciascun 
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coefficiente h la somma dei coefficienti omologhi delle prime due maMd 
nienzionate. In diversa forma di aWone, si pud prevedere di utto 

stimaton locali, i prodotti degli stimaton locaH raccolti in due o piii matrici di 
stimaton locaH. Si puo ipolizzare anche I'ulilizzo di tutti i coefficieati del polinomio 
5 interpolante. Nel caso di un polinomio del quarto ordiue si polx^ in tal caso 
combinare tia loro (sommandoli, moltiplicandoU o sottopouendoU ad altro tipo di 
elaborazioue) cinque coeffidenti otteaendo. da dnque matrid una matrice singola. 

Tomando all'esempio illustrato, che usa un singolo coefficiente dd 
polinomio interpolante e quindi una singola matrice. partendo dalla matrice degU 
10 stimaton locaHdiFig.7F,vienedetenninatol'istograiiima 
degU stimaton locali stessi in varie classi di appartenenza 

Sulla matrice degli stimaton locali considerata viene fetta scorrere una 
. finestra di dim^oni I x J opportunamente definite. La finestra viene posizionata in 
una prima posizione (Fig.7G) e in questa posizione vidie eseguito il calcolo 
15 deil'istograimna di distribuzione dd coeffidenti ay della matrice degli stimaton 

locali che ricadono nella finestra. 

Come si h osservato con riferimento alia Fig. 10 e seguenti, relative ai risultali 
sperimentaH, gli istogrammi relativi alle distribuzioni dd coeffidenti dd van ordini 
di un dato polinomio variano assumendo fonne diverse a seconda del tipo di tessuto" 
20 che si trovanella specifica poizione di immagine esaminata, raccMusa dalla fineslxa. 
Una volta posirionata la finestra di indagine di dimensioni IxJ in una data posizione 
sulla malrice de^ stimaton locali, si esegue I'istogramma di distnWone e si 
calcola un fettore di fonna dell'istogramma stesso, ad esempio la "deviazione 
standard" a. Spostando di un passo per volta la fikestra sia lungo le righe che lungo 
25 le colonne della matrice degU stimaton locali. si ottiene una serie di valori di 
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deviazioBe staaadard pari la nmnero di coeffidfeoti della malrice degli stiniatori locaH: 
mdicando con ag U valoxe deUa deviadone standard relativa aU'istogramma di 
distriburione dd coefficienti della matrice degU stiinatori locaH cdcolato con la 
fexestra in posizione ij. e possibile costruire una matrice con m nxnnero di 
5 coeffidentiparialnimeroddcoefficien1ia8dellan.aMcedegUstimatoril^^ 

cui valori sono stati modificati ed elaborati Imeado conto del parametro di fonna 
(nell'esetnpio spedfico la "deviazione standard" dell'istogranmm di 

distribuzione. Tale matrice sarii costituita da mxn valori bg 

bg = fl^ag, ag) 

10 dovef^unafonzionequalsiasidiay e ag. Ad esempio, si-pu6prevedere di divideie 
dascun coeffidente ag della matrice di stimatori locaU per il corrispondente valore di 
dij ottenendo: 

ba = aij/aij 

Si deve. peraltro, compx^dere che il tipo di funzione che lega i coefficienti by ai 
15 coefficienti ag ed al fettore di fonna ag dell'istx>gramma di distribuzione 
corrispondente pn5 essere diversa da.quella indicata. Ad esempio si pu6 prevedere di 
eseguire una moMpUcazione anzidi6 una divisions, oppure una diversa 

combinazione di qiiesti valori. 

Ancora, si pu6 ipolizzare di considerare per una plural di coefficieati di 
20 oxdine diverso del polinomiointerpolante una corrispondente 

stimatori locali e di combinare questi tea loro pesandoU con i parametri di 
lelativiistogrammidi distribuzione. |^ 

Anche il fettore o paramelro di forma dell'istogramma che viene\^li ^1o ^ 
pn5 essere diverse dalla. deviazione standard. Ad esempio si puo l^^^p^ 
25 utilizzare la curtosi, I'indice di simmetria, momenti di gmdo ^^f^^ 
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parainelro indicativo delie classi piil popolate. Ib geaerale, pu5 essere prevkto di 
elaborare i coeffideoli deUa matrice degli stimatori locaU con m qoialmxque 
paxaxnetrodifonnachecaratlmzzamogranmiadidi^^ . 
tale xnaMceia urn detenninata finest^ 
5 singolo coefficiente dellamatrice degli stimatori locali. 

La matrice di dimensione m x n che laccoglie i coefficienti h, viene 
dmominata nel seguito malrice degli stimatori locaH pesali. Ci5 sta a significare chc 
ciascmi coefficiente di tale matrice stato "pesato" per nn fettore di forma 
dell'istogramma di distribuzione della finestra corrispondente, dove per "pesato" si 
10 deveintendereilMtochetalecoefficierrte^statoinq^chemodocombinato 
parametro di forma, come la deviazione standard od altro parametro idoneo. 

Si deve comprendere che I'operazione descritta pu6 essere eseguita 
antonomamentc per tutti o per alcuni dei coeffidettti dci vari ordini del polinomio 
• interpolante che risultano utiH per la caratterizza^ione spetlxale del segnale 
15 «x,grafico. L'operarione eseguita sulla singola matrice degU stimatori locali che 
raccogUe i coeffidenti di m generico ordine del polinomio inteipolante pu6 essere la 
stessa per tutti gli ordini considerati, ma non si esclude di elaborare in modo diverso i 
coeffidenti dei vari ordini dd polinomio inteipolante. Ad esempio, per i coefficient! 
ddterzo ordine pu6 essere previstannapesatuialramite divisione per il paramdxo di 
20 deviazione standard, mentre per i coeffidenti dd secondo ordine pu5 essere prevista 
la pe^ od altm operazione di comWione con nn parametro diverso. quale la 
curtosi, oppure nna operadone diversa dalla divisione tm stimatore locale e 

parametro di forma, 

n concetto generale 4, comnnque, quello di combinare il vaiore dd 
25 coeffidenti rappresentanti lo stimatore locale per un coefficiente che tiene conto 
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■ deUa fonna deU'isto^amma di distribuzione dei coefficieati stesd inima determmata 
finestra, 

la Fig.7J h riportata schematicamente la procedura clie pu5 essere eseguita 
quando zrmcb6 m singolo coefficieate del polinomio interpolante vengono utilizzati 

5 tutti i coeffidenti di tale polinonuo per ottenere m ^ complesso ed elaborato 
insieme di stimatori locaU pesati. La matrice degU stimatoii locaU h quindi oma 
matrice liidimensionale di dimensioni m * n * k, dove m 6 il nmnero di campioid per 
ixaccia, n 6 il nmnero di Hacce del ftame ecografico digitalizzato e campionato e k b 
il numero di coefficiend del polinomio interpolante (polinonuo di grado k-1). Con 

10 ai/i\...aij<'^sono indicati i coeffidenti della matrice. . 

•I coeffidenti 4'\ -aif > vengono tra loro combinati per ottenere nn 
coeffidente c« di una matrice composta di stimatori localL Questa matrice 6 una 

matrice bidimensionaledi dimensioni n*m. Ciascuncoeffidentediessaddatada 

Cij = W>.aa<^\ ...a«*^ 
15 con f&nzionegenerica. Ad esempio pud essere una sommatoriaodunprodottodi 

tutti i k coeffidenti di posizione ij dellamatrice tridimensionale. 

Su questa matrice cg si calcolano i coeffidenti di foima a«, ad esempio la 
dfiviazione standard, degli istogrammi di distribuzione 1ramit« la procedura gi^ 
descritta con riferimento alle Figg.7G, 7R Dai coeffidenti Cij e dai parametri di 
20 forma ag si ottiene nd modo gi^ descritto la matrice degU stimatori locaH pesati, 
ancora indicati con bg e da questa la matrice di popolaaone fg. 

La procedura h concettualmente simile a quella descritta con riferimento alle 
Figg.7D-7H. con la differenza, tuttavia, die nella matrice finale sono contenuti 
coefficienti (stimatori locaU pesati) che contengono informazioni derivanti da tutti i 
25 coeffidenti dd polinomio interpolante. anzidi6 da uno solo di essi. 
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Sono anche possibili.soluziom intennedie, dove si utilizzano solo alcuni dei 
. coe£5cieati del polinotoio interpolante. 

In una variante di lealizzazione, si pu5 prevedere di eseguire prima la 
pesatura dei singoK coeflBcienti atteaverso il coefficiente di forma CTg e di eseguire 
5 solo successivamente la combinazione degU stimatori locali pesati tra loro. In 
sostanza, in tal caso detemiinata la matrice dei coefficienti spettraU si detemiinano i 
polinomi inteipolanti e dai coefficienti di questi la matrice degli stimatori ay*\ Per 
ogni valore di k si determina al distribuzione statistica, cioh la forma degU 
istogrammi e se ne calcola un fettore di forma, ad esempio la deviazione standard 
1 0 CTy®. Successivamente si calcolano gli stimatori locali pesati 

bij<^> = ^ay*>.<T/>) 
Infine, si combinano gli stimatori locali pesati per ottenere 

Cij = f(b8« ....bg*)) 
dove f e f sono fimzioni qualsiasi. 

1 5 Secondo un diverso aspetto dell'invenzione, il procedimento non prevede la 

pesatura dei coefficienti o stimatori locaK. Questa diversa forma di attuazione del 
procedimento non consente di utiUzzare le informazioni contenute nella fonna degli 
istogrammi di distribuzione, ma consente ancora di utilizzare il contenuto 
informative di pit coefficienti del polinomio interpolante. La Kg.7K rappresenta in 
:0 modo schematico la procedura di elaborazione in questo caso. SimboK uguali 
indicano elementi uguaU od equiyalenti a quelU deUa Fig.7J. La matrice c^ non viene 
piii elaborata per ricavare il fettore di forma cTij, bensi i suoi coefficienti (stimatori 
locali, non pesati) vengono sottoposti ad analisi statistica per ricavare direttauiente la 
matrice di popolazione, qui indicata con fij*. 

Ottenuta la matrice degU stimatori locaH pesati, i valori dei coefficienti di tale 
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matrice veagono classificati secondo classi di valori che sono risultate, 
speiirQentalmente, legate in modo biimivoco ad ima detenninata strattura tissutale. la 
altri termim\ si verifica sperimentalinente che m certo tipo di tessuto (ad esempio il 
tessuto neoplastico della prostata) altera 11 segnale ecografico di eccitazioiie in naodo 
5 tale che m determinato stimatore locale pesato assume valori che ricadono in una 
detenninata classe di valori. L'ecografo viene quindi programmato in modo tale da 
eseguire, su un frame in ingresso, le operazioni sopra descritte fino all'ottenimeato 
della matrice degli stimatori locali pesati. Quindi, i coelBficienti della matrice cosi 
ottenuta vengono confrontati co^ la classe di valori caratteristica del tessuto 

10 neoplastico e suU'immagme B-mode in bianco e nero viene sovrapposta una 
immagine cromatica formata dai punti della matrice degli stimatori locali pesati i cui 
valori ricadono nella classe caratteristica del tessuto neoplastico. Questa indicazione 
cromatica fomisce all'ecografista ririformazione richiesta, mettendo in evidenza la 
presenza e la'posizione del tessuto neoplastico. - ^ 

15 Si comprendera che la classificazione pud essere fatta per vari tipi di tessuti o 

di strutture. Ad esempio, se si desidera evidenziare la presenza di tessuto sano, sara 
sufficiente individuare la classe di valori di un certo stimatore locale pesato che e 
caratterisrtica della presenza di tessuto sano. Analogamente, quando si desideri 
evidenziare sulla stessa immagine simultaneamente o sequenzialmente la presenza di 
20 due diversi tipi di tessuto o di strattura (ad esempio la placca lipidica e la 
calcificazione nella carotide), b sufScieute individuare le due classi di valori che 
caratterizzano questi due tipi di stratture e quindi filtrare i coefScienti della matrice 
de^ stimatori locali con i valori caratteristici di queste classi e presentare il risultato 
simultaneamente od in sequenza suU'immagine ecogra£ca B-mode, eventua 
25 preferibilmente con due colori diversi. 
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Per individuare sperimentalmente le classi di valori che pear un certo stimatore 
locale caratteiizzano un determinato tipo di tessuto si puo procedere esaminando ua 
campione sperimentale in cui, tramite analisi istologica, si individuano le strutture di 
interesse, Nel caso della prostata, ad esempio, si individuano le zone omogenee di 
6 tessuto sano» di adenocaidnoma e di adenofibroma. Quindi, si esegue relaborazione 
del segnale ecografico in radiofirequenza, ottenuto in risposta daUa scansione del 
campione sperimentale, come sopra descritto, fino ad ottenere una singola matrice 
degli stimatori locali pesati. 

Per come la matrice degli stimatori locali e stata ottenuta, ad ogm coefficiente 
1 0 della matrice corrisponde un punto dell'immagLne B-mode.'Perlanto, essendo noto in 
quali zone si localizza una determioata struttura, ad esempio radenocarcinonaa, jsi 
Verifica la distribuzione statistica degli stimatori locali in tali zone e si individua la 
classe con la maggiore firequenza di popolazione. Quella classe verra associata, per 
quel datb oigano alia presenza di adenocarcinoma. Pertanto, ogni qualvolta venA* 
15 eseguita una analisi ecografica su un organo analogo, gU stimatori locali pesati 
verranno classificati e quelli che ricadranno neUa classe di riferimento 
dell'adenocarcinoma verranno evidenziati sovrapponendo aU'immagine B-mode una 
scala cromatica dove risultano colorati i pixel corrispondenti alle zone dell'inunagine 
in cui i valori degli stimatori locali pesati ricadono nella classe caratteristica 
20 dell'adenocarcinoma. Analogamente si proceder^ per gli altri tipi di tessuto (tessuto 
sano, adenoiahroma, nel caso della prostata). 

Una elaborazione pit. articolata prevede di lavorare su almeno una coppia di 
roatrici di stimatori locali sia in fese di individuazione delle classi di appartenenza, 
sia in fase di analisi di un organo sottopo^to ad indagine ecografica. Un esempio di 
!5 elaborazione di questo tipo, basata su tre stimatori locali distinti, viene oira descritta 
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con r ausilio deUe Figg.23A e 23B. Sul campione sperimentale viene ancora es^ta 
I'acquisizione di rm fiame ecografico, che viene poi elaborato come sopra descritto 
fino aU'ottenimento di tre matrici di sftimatori locali pesati. 

Su nn primo grafico bidiinensionale si riportano, in asdsse ed in ordinate 
5 rispettivamente, tutti i valori dei coeffidenti di due matrici di stimatori locaU pesati, 
ad esempio quelli ottenuti dai coefficienti di ordine 0 e di ordine 1 del polinomio 
inteipolante. Su un secondo grafico si riportano gU stimatori locali pesati ottenuti dai 
coefficienti di ordine 0 e di oidine 2. In Fig.23 A sono riportati in ascisse i coefficienti 
della matrice degU stimatori locaU pesati otteunti dai coefficienti del primo ordine del 
1 0 polinomio inteipolanfe per una prostata. In ordinate sono riportkti i coeffidenti della 
matrice degU stimatori locaK pesati ottenuti dai coeffidenti di ordine 0 dello stesso 
polinomio. In Kg.23B sono riportati in ascisse gU stimatori locali pesati ottenuti dai 
coeffidenti del secondo ordine ed in ordinate queUi ottenuti dai coeffidenti di ordine 
.0. 

15 Ad ogni punto di dascuno dej due grafici risulta co^ associata una coppia di 

valori Oe coordinate del punto) costituite dai valoie" deU'uno e dell'altro degU 
stimatori locali pesati utilizzati. 

I grafici die si ottengono in questo modo sono sovrapponibiK al j&ame 
ecografico in ingresso. Si possono cosi associare ad ogni punto dd fiame ecografico 
20 di ingresso coppie di classi di valori degU stimatori locali pesati piesi in esame e 
conispondenti a tale posizione. 

Sul fi:ame ecografico in ingresso si individuano, tramite analisi istologica, le 
zone omogenee per natura dd tessuto (sano, neoplastico etc.). Nelle Figg.23 A e 23B 
queste zone sono contraddistinte da tre tonality di grigio diverse. Poich6, in entrambi 
55 i grafici, ad ogni punto del fiame ecografico corrisponde una coppia di valori degli 
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stimatori locaU pesati, divieae possibile associare ad ogni zona omogenea del tessuto 
esaminato una classe di valori per detti stimatori locaK pesati, classe univocamente 
Msociata a quel detemiinato tipo di tessuto. 

Per poter individuare biunivocamente ogni zona 6 importante selezionare un 
5 insieme costituito daUe coppie o triplette di classi di valori degli stimatori locali 
pesati. taH da non sovrapporsi con nessun altro associate ad una differente zona. 
L'operazione di estrazione di taH insiemi pud essere fetta attraverso operazioni 
logiche di tipo OR e AND, ad esempio. 

Come accennato, la fiequ^iza di popolazione degU insiemi cosi individuati 
10 vienerappresentataconunascalacromatica,diversaperogbiiiisieme,suU'imm^ 
B-mode tradizionale ed b direttamente proporzionale aU'intensiti della fiequaiza di 
popolazione. 

Sperimentalmente, come peraltro si pud notare anche dalle Figg. 17A, 17B, 
normalmente I'aiea che interessa evidenziare (nell'esempio il tessuto neoplastico) d 
15 evidenziabile combinando per sempUce somma le immagini ottenute da tutti i 
coefficient! del polinomib intetpolante, evidenziabile combinando, in questo caso, 
con un'operazione di AND gU insiemi di valori ottenuti da tutti i coefficienti del 
polinomio intetpolante. 

La proceduia sopra descritta per la correlazione sperimentale tm 
20 caratteristiche tissutaK rilevate tramite analisi istolbgica e classi di valori per gli 
stimatori locali pesati vale anche nel caso di stimatori locali non pesati, ottenuti ad 
esempio combinando tra loro coefficienti di diverso gtado dei polinomi iuteipolanti. 

Nei lisultati sperimentali sopra discussi ed esposti nelle Figg.9A-17B si h 
utilizzata una decomposizione del segnale RF in sedici bande o sottobande di uguale 
55 larghezza. Questa scomposizione ha dato risultati utiH tanto neU'esame della 
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prostata, quanto neU'esame deUa carotide. Questi due tipi di organo sono stati sceW 
in quanto presentano strattuie tissutaH di dimensioni diverse. La ghiandola piostatica 
presenta tipicamente struttuie di dimensiom maggioii lispetto alia caiotide. La 
diversa dimeasione deUe strutture che riflettono U segnale ecografico di eccitazione 
5 comporta una diversa distribuzione delle fiequenze nel segnale di ritomo. La 
• suddivisione in sedici sottobande di uguale ampiezza consente di ottenere una 
lisoluzione sufficientemente significativa del segnale RP restituito da entrambi gU 
organi esatninati. 

In generaie, tuttavia, si comprender^ da quanto sopra esposto che la- scelta 
10 della decomposizione in bande (numero di bande, loiii ampie2^) pu5 essere diversa 
a seconda del tipo di tessuto o di struttura da analizzare. Tanto pit piccole sono le 
strutture che si desidem studiare neU'oigano analizzato, tanto maggiore dovr^ essere 
la risoluzione della decomposizione neUe alte fiequenze e viceversa. Ci6 h 
confennato dai risultati sperimentaK di seguito discussi con riferimento aUe Figg. 1 8 a 
15 21, ottenuti ancora dall'elaborazione di una immagine ecografica di una porzione di 
carotide. 

In Figg.l9, 20 e 21 6 riportata I'immagine B-mode della porzione di carotide 
analizzata. fl segnale in RF, digitalizzato e campionato, h stkto decodificato tramite 
una wavelet discreta in sei bande di ampiezza variabile. In Fig.l8A d mostrato 

20 I'albeio r^ppresentativo del process© di filtraggio seguito ed in Fig.lSB le sei bande 
di decomposizione. Si osserva che I'ampiezza deUe bande cresce daUe fiequenze 
basse alle fiequenze alte. Con fc h indicata la fiequenza di campionameuto del 
segnale. In Fig.l8C d mostrato I'andamento del polinomio interpolante ottenuto '-^ 
applicando la decomposizione suddetta all'immagine B-mode delle Figg. 

*'5 In Fig.l9 6 riportato: nella parte superiore I'immagine B- mode 4^U?^^otide 
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esaminata; al di sotto di essa I'andameato del polinomio interpolante del quarto 
ordine (comspondente alia Fig.l8C) e successivamente gU istogtammi di 
distribuzione dei coefficienti A0-A4 dei van gradi del polinomio intopolante, 
calcolati in corrispondenza di una finestra (mostrata in sovrapposizione all'iimnagine 
5 B-mode) centrata suUa placca lipidica. In Fig.20 con analoga disposizione delle 
immagini sono riportati il polinomio inteipolante e gU istogrammi di distribuzione 
ottenuti in una finestra centrata suUa zona Ubera del vaso, ciofe sul sangue. Anche in 
questo caso la finestra h mostrata in sovi^>posizione aU'immagine B-mode della 
carotide. Infine, in Fig.21 h mostrata, con analoga disposizione, U risultato ottenuto in 
1 0 una finestra centrata su una porzione.dell'immagine B:mod-e in bui si xiconoscono dei 
fenomeni di calcificazione. 

Dal confionto deUe Figg.l9, 20 e 21 si osserva che con la decomposizione in 
sei bande di ampiezza crescente daUe basse aUe alte j&equenze I'andamento del 
polinomio inteipolante non cambia in nibdo apprezzabUe nelle tre zone osservate 
15 Cplacca Upidica, sangue e calcificazioni). Questo b un piimo indice del fetto che la 
decomposizione adottata non fomisce una lisoluzione idonea per I'esame di questo 
tipo di oigano. Ci6 b dovuto al fetto che le strutture che devono essere riconosoiute 
sono di piccole dimensioni e quindi I'infoimazione spettrale chte esse fomiscono si 
ritrovano alle alte firequenze dello spettro di risposta, ciod Ik dove la decomposizione 
20 adottata prevede bande di iampiezza maggiore. Contrariamente a ci6, i tisultati 
ottenuti suUo stesso organo con una decomppsizione in sedici bande di uguale 
ampiezza (Figg.9A a 16) moslrano che la lisoluzione ottenuta con tale 
decomposizione h sufficiente a fomire risultati significativi. 

DagK istogrammi di distribuzione deUe Figg.l9, 20 e 21 si nota che in effetti 
!5 - bench6 I'andameato del polinomio inteipolante sia sostanziahnente uguale neUe Ire 

Pag.38 di 44 

C:\Documeitfs and ScttingsVinichete tnannucciCesktopWocuinenti inichele\actis_081003\iict300903def (2).doc 



zone di placca Kpidica, sangue e calcificazione - la foima degU istogrammi di 
distribiizione dei coefficienti dei vari ordini cambia. Tuttavia, la differ^nza tra 
istogrammi di distribuzione di coeffidenti omologhi nelle tre zone indagate non h 
altrettanto marcata e significativa di quella che emerge utilizzando ima elaborazione. 
6 del segnale basato su decodifica in sedici bande di uguale ampiezza. Ih 

definitiva, il confionto tra i risultati ottenuti con i due tipi di decodifica in sottobande 
dimostra che la scelta di una decodifica opportuna, per numeax) di bande e ampiezza 
delle stesse. consente di adattaxe il metodo di anaHsi della presente invenzione per 
ottenere risultati ottimaU a seconda del tipo di tessuto o di struttura indagata. La 
10 scelta del tipo di decodifica o scomposizione pud aweniie sperimentalmente, cosi 
come sperimentalmente si deteiminano il grado del polinomio intetpolante piii 
adatto, i coefficienti di esso piCi significativi ed i coefficienti di fonna piii lilevanti 
degli istogrammi di distribuzione di detti coefficienti. 

In Fig. 22 h mostrato esenq)lificativamente un albero di decomposizione dello 
15 spettro del segnale in radiofrequenza tramite la Wavelet discreta DWPT. A seconda 
del tipo di indagine ecografica che deve essere eseguita e delle caiatteristiche della 
struttura esaminata, il metodo di indagine pud utilizzare una qualunque delle vaiie 
possibility di decomposizione lungo Talbero di Fig.22. 

E' inteso che il disegno non mostra che una pratica forma di attuazione 
20 dell'invenzione, la quale pud variare neUe forme e disposizioni senza petaltro uscire 
dall'ambito del concetto alia base dell'invenzione. L'eventuale presenza di numeri di 
liferimeato nelle accluse rivendicazioni ha il solo scopo di feciUtame la lettura alia 
luce del teste che precede e degli aUegati disegni e non ne limita I'ambito di 
protezione. 

25 
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RIVENDICA7TQNT 
1. Un metodo di analisi spettrale di un segnale ultrasonico a 
ladiofiequenza restituito da una struttura sottoposta ad un esame ecografico, 
comprendente le fiisi di: 

5 a) trasmettere un segnale ultrasonico di eccitazione in una porzione di detta 
struttura sottoposta ad esame; 

b) ricevereun segnale a radioJ&equenzadirisposta da detta struttura; 

c) appHcare una sequenza di filtraggi in modo da ottenere una decomposizione 
della banda del segnale a radioficequenza di risposta in una pluraliti di bande di 

10 fiequenza; "... 

d) dai coefficient! risultanti da detto mtraggio calcolare stimatori locaU (ag, cg), 
contenenti informazioni sullo sfpettro del segnale a radiofrequenza; 

caratterizzato dal fatto di c ombinare detti srtunatoii locaU (ay; cg) con parametri 

(a„) rappiesentativi della fonna di una distribuzione statistica di detti stimatori locali 

15 in una porzione di una immagine ecografica. 

2. Metodo come da rivendicazione 1, in cui le bande di fiequenza in cui 
viene suddiviso detto segnale a ladiofiequenza sono taH da ricoprire I'inteia banda di 
fiequenza del s^nale. 

3. Metodo come da rivendicazione 1 o 2, in cui dette bande di fiequenza 
20 sono bande di diversa larghezza e posizione. 

4. Metodo come da una o piu deUe livendicazioni precedent!, 
comprendente le fasi di: 

- per un fiame ecografico in ingresso generare un fiame campionato e 
digitalizzato; 

15 - decomporre detto fiame campionato e digitalizzato in dette bande di fiequenza; 
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- generare una matrice di coeffidenti spettrali contenente i coefficienti lisultanti 
da detto jSltraggio o coefficienti da essi derivanti; 

- determinare, per alcuni almeno dei pimti del fiame campionato e digitalizzato, 

rispettivi polinomi inteipolanti (PI) che approssimano I'andamento di detd 

5 coefficienti spettraU neUe varie bande in cui 6 state decomposto H segnale a 

radiofrequenza; 

- per detti punti, da almeno tmo dei coefficienti (ao,...ak) del polinomio 
interpolante ricavare detto stimatore locale (ajj; c^), detti stimatori locaU 
costituendo una matrice di stimatori locali. 

^- Metodo come da rivendicazione 4, in cui ciascuno di detti stimatori 
locaU (aa) h costituito da uno dei coefficienti del rispettivo polinomio interpolante. 

6. Metodo come da rivendicazione 4, in cui per dascun punto vengono 
detenninati almeno due stimatori locaU (aQ^>) suUa base di almeno due coefficienti 
del polinomio interpolante, per generare una matrice tridimensionale di stimatori 

15 locali (ag^J) . 

7. Metodo come da rivendicazione 4, in cui ciascuno di detti stimatori 
locali (cij) h costituito da una combinazione di una plurality di coefficienti del 
cortispondente polinomio interpolaiite. 

8. Metodo come da rivendicazione 4 o 5 o 6 o 7, in cui ciascuno di detti 
20 stimatori locaU (ag; cy) viene combinato ad un coefficiente di fonna (cij) di un 

istogramma di distribuzione di detti stimatori locaK in una finestra entto cui d 



IMARCADaBOLLO 



contenuto detto stimatore locale, per ottenere uno stimatore locale pesato (bijy^^^i 
9. Metodo come da rivendicazione 8, comprendente le fasi di:/y iS> '. 
• determinare una disfaibuzione statistica di detti stimatori locali •) in J 

25 finestre di dimensioni inferiori rispetto alia dimensione di detta matricB^l. ., 

Pag.41 di44. IE 

CiVDocumenls and SettiiigsVmichele niaiuiucci\Deskipp\docuinenti iiiichde\8ctis_0810Q3\act30Q903def ^) doo 




stimatori locali (^j, C|j); 

• detemunare un parametro di forma (ag) caratteristico di detta distribuzione 
statistica per ciascuna di dette fuaestre; 

• per ciascuna finestra, combinare detto parametro di fonna (CTg) con m 
5 conispondente stimatore locale (ag;ca) per ottenere imo stimatore locale pesato 

(bij). 

10. Metodo come da una o pitl delle rivendicazioni 4 a 9, in cui piii 
stimatori locali pesati ottenuti per uno stesso punto del fiame ecografico campionato 
e digitalizzato tramile diversi coefficienti (ao, ... a^ del rispettivo polinomio 

10 intetpolane vengono tra loro combinati. 

11. Metodo come da una o piii delle rivendicazioni precedent!, in cui detti 
filtraggi sono ottenuti tramite tma trasformata tempo-fiequenza. 

12. Metodo come da rivendicazione 11, in cui detta trasformata tempo- 
frequenza d una Wavelet 

''^ . Metodo come da rivendicazione 11, in cui detta liasfonnata tempo- 

frequenza 6 una trasformata Wavelet Packet discreta (DWPT). 

14. Metodo come da una o pit deUe rivendicazioni precedenti, 
comprendente la fese di determinare la distribuzione statistica degli stimatori locaK 
pesati e di creare un insieme di classi di valori in grado di individuare 

20 biunivocamente porzioni omogenee sul frame ecografico del campione investigato. 

15. Metodo come , da una o pit delle rivendicazioni precedenti, in cui su 
una immagine ecografica vengono sovr^poste immagni cromatiche generate tramite 
detti stimatori locali pesati. 

16. Metodo come da rivendicazione 15, in cui dette immagini cromaticlie 
25 vengono generate selezionando gH stimatori locaU pesati che ricadono in classi di 
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appartenenza, biumvocamente legate a strutture tissutaK predetemmiate. 

•* 

I 17. Un metodo di analisi spettrale di un segnale tiltrasonico a 

radiofrequenza restituito da tma stnittuta sottoposta ad un esame ecografico, 
comprendente le fasi di: 
5 a) trasmett^ un segnale ultrasonico di eccitazione in una porzione di detta 
struttura sottoposta ad esame; 

b) riceveie in ingresso un segnale a r adiofrequenza di risposta da detta struttuia; 

c) per un frame ecografico in ingresso generate un frame campionato e 
digitalizzato; 

10 d) applicare a detto frame campionato e digitalizzato una sequenza di jBltraggi in 
modo da ottenere una decomposizione della banda del segnale a radiofiequenza 
di risposta in pna pluralita di bande di frequenza; 
e) generaie una matcice di coefficienti spettrali contenente i coefficienti risultanti 
da detto filtraggio o coefficienti derivanti da essi; 
15 f) determinare, per alcuni almeno dei punti del fi:ame campionato e digitalizzato, 
rispettivi polinomi interpolanti (PI) che approssimano Fandamento di detti 
coefficienti spettrali nelle varie bande in cui e state decomposto il segnale a 
radiofrequenza; 

g) per detti punti, dai coefficienti (ao,...^ del polinomio interpolante ricavare 

20 uno stimatore locale (cij) combinando tra loro almeno due coefficienti di 

diverse ordine (ao, ...Bk) d^l polinomio interpolante. 

18. Metodo come da rivendicazione 17, in cui le bande di frequenza in cui 

viene suddiviso detto segnale a radiofiequenza sono tali da ricoprire Tintera banda di 

< 

frequenza del segnale. 

25 19. Metodo come da rivendicazione 17 o 18, in cui dette bande di 
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firequaoza sono bande di diversa larghezza e posizione. 

20. Metodo come da una o piii delle livendicazioni 17 a 19, in cm detti 
filtraggi sono ottenuti tramite una trasformata tempo-fiequenza. 

21. Metodo come da rivendicazione 20, in cui detta trasfonnata tempo- 
5 fiequenza h una Wavelet 

22. Metodo come da rivendicazione 20, in cui detta trasformata tempo- 
firequenza 6 una trasformata Wavelet Packet discreta (OWPT). 

23. Metodo come da una o piii delle rivendicazioni 17 a 22, comprendente 
la fase di determinare la distribuzione statistica degli stimatoii locali pesati e di 

1 0 creare un insieme di classi di valori in grado di individuare biunivocamente porzioni 
omogenee sul jQrame ecograSco del campione investigato. 

24. Metodo come da una o piii delle rivendicazioni 17 a 23, in cui su una 
immagine ecografica vengono sovrapposte immagini cromatiche generate tramite 
detti stimatoii locali pesati. 

^ 5 25. Metodo come da rivendicazione 24, in cui dette iTmnag iTii cromatiche 

vengono generate selezionando gli stimatori locali che ricadono in classi di 
^artenenza, bhmivocamente legate a stnttture tissutali predetermlnate. 

26. Un dispositivo ecografico comprendente una sonda ecografica, mezzi 
di acquisizione e di elaboiazione un segnale di ritomo a radiofrequenza da una 
20 struttura sottoposta ad un esame ecografico, caratteiizzato dal fetto che detti mezzi di 
acquisizione ed elaixoazione sono programmati per eseguire un metodo come da una 
o pia delle rivendicazioni 
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